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The local quasiparticle density of states around a unitary impurity with a Kondo-like magnetic 
moment induced at its nearest neighbors in a d-wave superconductor is studied within the slave- 
boson mean-field approach. The Hamiltonian is exactly diagonalized and the results are obtained 
self-consistently. We show that the interference between the strong impurity potential scattering 
and the Kondo effect leads to novel quasiparticle spectra around the impurity, which are strikingly 
different from the case of a single unitary or magnetic impurity. The recent STM image of the local 
differential tunneling conductance around a Zn impurity in a high-Tc cuprate can be essentially 
explained if the blocking effect of BiO surface layer between the tip and probed Cu02 plane [Zhu 
et al, Phys. Rev. B 62, 6207 (2000)] is taken into account. 
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The scanning tunneling microscopy (STM) turns out 
to be a powerful probe of the local effects of individual im- 
purity atoms or defects on the fundamental properties of 
high-Tc cuprates [0H| . The experiments revealed the in- 
duction of a virtual bound quasiparticle state around the 
individual impurity, providing a strong evidence of the 
d-wave pairing symmetry of high-Tc superconductors . 
This observation is motivated by the earlier predictions 
of Balatsky, Salkola and co-workers H,^. However, the 
high resolution STM imaging of the effect of individual 
impurity atoms Zn substituted with a controlled man- 
ner for Cu in BizSrzCaCuaOg+a (BSCCO) [| also pre- 
sented an unexpected spatial distribution of the local 
density of states (LDOS) strength at the resonance en- 
ergy: It has a strongest intensity directly on or above 
the Zn site, local minimum on its nearest-neighboring 
sites and local maximum on second-nearest-neighboring 
sites. This distribution is opposite to all theories based on 
the strong atomic-like nonmagnetic impurity model , 
which show a vanishing density of states directly at the 
impurity site and a strong resonance peak on its near- 
est neighbors. Recent NMR measurements of Zn and Li 
doped YBa2Cu306+j, |p|-pd]| and the specific heat mea- 
surement of Zn-doped YBa2Cu3 6.95 indicated that 
a Kondo like magnetic moment with S* = ^ is induced at 
the nearest neighboring sites of the Zn (or Li) impurity. 
Motivated by the latter experiments, the Kondo effect on 
the local quasiparticle states around an induced magnetic 
moment in a d-wave superconductor has been theoreti- 
cally studied by several groups p3|-p^ . All these results 
show wide-distanced double-peak structure with respect 
to the zero bias at the impurity sites, and therefore do 
not agree with the STM imaging of the effects of a Zn 
impurity. We believe that these studies seem to be more 
suitable to describe the Kondo impurities such as Mn 
and Co substituted for Cu in high-Tc cuprates instead 
of Zn or even Ni whose spin is believed to have a fer- 
romagnetic coupling with conduction electrons [^-p^. 
In order to compare with experiments, the scattering of 



quasiparticles due to both Zn and the induced magnetic 
moment around it needs to be carefully investigated. So 
far this problem has not been properly addressed, par- 
ticularly in view of the scattering strength of a Zn im- 
purity being in the unitary limit as it will be estimated 
below. Valence counting suggests that the Zn impurity 
maintains a nominal Cu^+ charge, the Zn^+ ["id}^] with 
5 = configuration indicates it would act as a nonmag- 
netic impurity. According to Ref. jl^, the second and 
third ionization energies for the Zn element are respec- 
tively 18 eV and 39.8 eV, while the third ionization en- 
ergy of the Cu element is 36.9 eV. In the cuprates, the 
relevant Zd Cu^"*" electrons form a narrow d^^-y^ band 
with bandwidth less than 2 eV. The center of the band, 
where the chemical potential lies, is located around - 
36.9 eV below the vacuum. The top of 3^^*^ level of 
a Zn^+ ion is at —39.8 eV which is below the bottom 
of the dx2-y2 band. The conduction electrons near the 
chemical potential will experience a repulsive local po- 
tential Vo « (-18 eV) - (-36.9 eV) = 18.9 eV at the Zn 
site. For such a large local potential as compared to the 
bandwidth, Zn can indeed be regarded as a unitary im- 
purity. Other evidence for the strong potential scattering 
from Zn has also been provided experimentally It 
is important to point out that our present consideration 
is very different from the case studied by the work p^ ] 
where the authors assumed the Zn impurity to have a 
weak attractive potential. 

In this Letter we study the local quasiparticle spectrum 
around a unitary impurity with a Kondo-like magnetic 
moment induced at its nearest neighbors in a d-wave su- 
perconductor. Based on this more sophisticated model, 
it is shown that the presence of the strong potential scat- 
tering from the nonmagnetic impurity leads to a stronger 
single resonant peak in the LDOS directly at the site of 
the magnetic moment, in contrast to the case for a sin- 
gle magnetic impurity [ld|Jl5[| . Comparing with the case 
of a single nonmagnetic unitary impurity, the LDOS is 
strongly enhanced by the Kondo effect from the induced 
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magnetic moment. We also show that the novel pattern 
observed in the STM experiments can be qualitatively 
explained by taking into account the blocking effect pl| ] 
of the BiO surface layer which exists between the tunnel- 
ing tip and the Cu02 plane where the Zn impurities are 
located. 

We model the d-wave superconductor defined on a 
two-dimensional lattice within a phenomenological BCS 
framework. The unitary impurity is taken to be at the 
origin Tq = (0,0). Nearest- neighboring to the impurity, 
the induced magnetic moment is described by the An- 
derson s-d exchange model. Since we are most inter- 
ested in the low temperature quasiparticle physics, the 
on-site repulsive interaction on the moment-sitting site 
is assumed to be infinite so that the double occupancy 
on this site is forbidden. We use the slave-boson mean- 
field approach [p2|j2^ ] to study the Kondo effect associ- 
ated with the induced magnetic moment. The system 
Hamiltonian can then be diagonalized by solving the the 
Bogoliubov-de Gennes equations [Q: 

with the bosonic number bo and the d-wave bond pairing 
amplitude Ay subject to the self-consistent condition 

Ay ^ ^ + ufvr )tanh{En/2kBT) , (2) 
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and 

bl = 1 -2j2{\u:ffiE„) + Kf[l fiE^)]} . (3) 

n 

f uV- \ . 

Here ( n 1 is the quasiparticle wavefunction with eigen- 
value En- Hij = -tij - flSi^nj + Kjf^iro + (erf + Ao)^iri(5ij 

is the single particle Hamiltonian. Without loss of gener- 
ality, we temporarily assume the magnetic moment to be 
at ri = (f , 0), S are the unit vectors of the square lattice, 
iij = Vbo for (ij) = (ri + S,ri) or (ri,ri + S) and 
t otherwise with V being the coupling strength between 
the moment and conduction electrons and t the hopping 
integral between conduction electrons, fi is the chemical 
potential, Vq is the strength of impurity potential, is 
the bare energy level of the moment electron with re- 
spect to the chemical potential for the band electrons for 
superconductivity, Aq is the Lagrange multiplier intro- 
duced to enforce the single occupancy constraint. Here 
f/ij represents the strength of d-wave pairing interaction, 
the value of which is gy — for (ij) = (ri + ^, ri) or 
(ri, ri + S) upon the assumption of different nature of 
the conduction and moment electrons and constant 
g otherwise. Finally, f{E) = [1 + ex^p^E/ksT)]-^ is the 
Fermi distribution function. 

Throughout the work, our numerical calculation is per- 
formed on a square lattice with size Nl = x Ny = 



35 X 35 and averaged over 36 wavevectors in the super- 
cell Brillouin zone. The strong scattering potential of 
the unitary impurity is taken to be Vb = 100, the bare 
energy level = —2, the chemical potential = —0.2, 
and the pairing interaction g = 1.2. Unless specified oth- 
erwise, hereafter we measure the energy in units of the 
hopping intergal t and the length of the lattice constant 
a. A trivial solution to Eq. (|^) with Aq = —erf and 6o = 
is always found, which describes the moment decoupled 
from the conduction electrons. However, the physically 
desired value of these two parameters are determined by 
minimizing the free energy with the procedure as detailed 
in Ref. ||lj|. With the chosen parameter values, we find 
a critical value of the coupling strength Vc « 0.9, be- 
low which the moment is decoupled from the conduc- 
tion band. Compared with the single Kondo impurity 
case, the larger Vc obtained here is due to the fact that 
the conduction electron density is strongly depressed on 
the unitary impurity site nearest-neighboring to the mo- 
ment so that the screening effect is weakened. When 
V > Vc, this free local spin state gives way to the Kondo 
screened state. For V = 1.0, it is found that Aq = 2.15 
and bo = 0.56. The self-consistent solution also gives 
a d-wave order parameter which is strongly suppressed 
around the unitary impurity in the scale of the coher- 
ence ^0 = TT'Vp /ir/S.g, where is the Fermi velocity and 
Ag « 0.4 is the bulk d-wave energy gap. 

Once the self-consistent solution is obtained, the LDOS 
is then calculated according to 

Pi ^ 2 ^[|<P^(£;„ -E) + \v-\H{E^ + E)] , (4) 

which is proportional to the local differential tunneling 
conductance as measured by the STM experiments. Our 
numerical result indicates that the LDOS is vanishing at 
the unitary impurity site, consistent with previous cal- 
culations by others 1^,0. Figure shows the LDOS as a 
function of energy on sites nearest-neighboring to the uni- 
tary impurity site, A = (1,0) [i.e., the magnetic moment- 
sitting site], B ^ (0,1), C = (-1,0), D = (0,-1) and 
on the corner site of the supercell. Due to the spatial 
inverse symmetry with respect to the co-linear axis con- 
necting the impurity and moment sites, the LDOS on 
sites B and D are identical. Notice that the LDOS on 
the corner site of the supercell recovers the bulk den- 
sity of states of a clean d-wave superconductor, which 
indicates that the supercell size and the number of wave 
vectors in the Brillouin zone has been large enough to un- 
cover the physics intrinsic to an isolated impurity with an 
induced moment. The LDOS on the nearest-neighboring 
sites has evidenced the existence of near-zero energy reso- 
nant states by showing the near-zero energy peaks. Com- 
pared with a single nonmagnetic unitary impurity or a 
single Kondo impurity, the unitary impurity with an in- 
duced moment shows much richer characteristics of the 
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quasiparticle resonant states in its vicinity: (i) In strik- 
ingly contrast to a single Kondo impurity case [ p^|jT^ ] , the 
LDOS on the magnetic moment-sitting site-A only dis- 
plays a strong single near-zero-energy peak, which comes 
in a unique manner from the presence of the unitary im- 
purity, (ii) Different from a single nonmagnetic unitary 
impurity case |^,||, the LDOS on site-B and site- 1? ex- 
hibits the splitting of the near-zero-energy peak. This 
difference is mainly due to the presence of the moment 
near the unitary impurity, which plays the role of weak 
scattering center and reduces the effective hopping inte- 
gral locally. As a comparison, we have also considered a 
(lO)-oriented dimer consisting of two nearest-neighboring 
unitary impurities and found that the splitting is much 
increased. Therefore, in the Kondo resonance regime, the 
unitary impurity plus a nearby moment composite can be 
effectively regarded a tailed nonmagnetic unitary impu- 
rity. Since the influence of the induced moment is very 
weak on site-C, the LDOS on this site only shows a sin- 
gle near-zero energy peak, (iii) As a result of the Kondo 
resonance, the intensity of the near-zero-energy peak on 
site-A is about four times large as that on other three 
nearest-neighboring sites of the unitary impurity. So far 
we have considered the local quasiparticle resonant states 
around the unitary impurity with a moment induced stat- 
ically. When a Zn impurity is substituted for Cu in 
a high-Tc cuprate, the 5* = ^ moment associated with 
Cu^"*" is removed from the system. As a result the local 
antiferromagnetic pairing is broken and a moment with 
spin-i is induced at any of the four neighboring sites of 
the impurity with equal probability. Our treatment here 
is consistent with the NMR experimental indication that 
the Zn impurity only reveals already- existing moments, 
localized on Cu sites of the doped antiferromagnet pC| ] and 
the induced spins on the nearest neighboring copper sites 
are not correlated psj . It is different from the approach 
based on a putatively single moment around the Zn im- 
purity like a spin cluster. After averaging over the four 
moment-position configurations around the impurity, the 
combined contribution from sites-A through D leads to 
the averaged LDOS p;"* = {pA + Pb + Pc + Pd)/^ at any 
of these four sites. As shown in Fig. ||, p?"* is character- 
ized by a very strong sharp peak close to the zero energy. 
Compared to the case of a single nonmagnetic unitary 
impurity, the corresponding peak intensity is enhanced 
by about 40%. Also interestingly, the LDOS spectrum 
also shows a narrow side-peak with a rather weak inten- 
sity. The highly asymmetric double-peak structure dis- 
tinguishes it with a sharp and narrow dip in between and 
their separation is only about 5% of 2Ag. In comparison, 
the double peaks obtained for a single Kondo impurity 
has a separation as large as 50% of 2 Ag . 

For the purpose to understand the LDOS spectrum 
and the spatial distribution of its strength at the res- 
onance bias observed in the STM experiment, we need 
to know the nature of the surface layer in BSCCO on 



which the tip is directly probing. There are ample ex- 
perimental evidences that the cleaved surface along the c- 
axis is always the BiO layer. The superconducting Cu02 
layer on which Zn impurities are located and the STM 
tip would like to explore, lies about 5A below the sur- 
face BiO layer. The crystal structure shows that each 
Bi atom sits vertically on the top of each Cu atom or 
Zn atom, and there is also an O atom associated with 
SrO layer between them. When the STM tip is trying to 
probe the Zn site, the hard core of the Bi'^"'" ion (with ra- 
dius 1.03 A) and also 0^~ (with radius 1.4 A) sitting on 
top of it will inevitably block the tunneling current 
from directly passing through the Zn site, where Zn^+ 
has the ionic radius 0.74 A. Because the STM image has 
such a high spatial resolution to the atomic scale, the 
tunneling current from the tip to the Cu02 plane is only 
distributed within a small area, the linear dimension of 
which is about one lattice constant. Therefore, the mea- 
sured LDOS above the probed Zn or Cu site is mainly 
contributed from its four nearest-neighboring sites, which 
can be approximately expressed as p°^^ = A^^p^^"^. 
Here A < 1 is a constant depending on the distance be- 
tween the tip and the four nearest-neighboring Cu atoms. 
As a result, the STM spectrum on the Zn site observed by 
experiment is in good agreement with p°^^ at i = (0, 0) 
which has the identical characteristic as that shown in 
Fig. ||. Experimentally, a small-side peak adjacent to the 
left of the strong resonant peak was indeed observed in 
the STM measurement ||l|. This is consistent with our 
result in Fig. 2 and it could very well be a signature 
of the induced Kondo-like magnetic moment around the 
impurity. In Fig. |3| a bubble plot is given for the spatial 
distribution of pf"^ at the near-zero-bias resonance energy 
E — 0.02, where the size of each black dot is proportional 
to the LDOS intensity at that lattice site. As it is shown, 
the pf"^ distribution agrees, at least qualitatively with the 
measured differential-tunneling conductance image ||l|]. 

To conclude, we have considered the effect of an indi- 
vidual nonmagnetic unitary impurity with induced mag- 
netic moment in a d-wave superconductor. To be applica- 
ble to the Zn impurity substituted for Cu in copper-oxide 
cuprates, the model has taken care of both the strong 
potential scattering from the nonmagnetic impurity it- 
self and the Kondo effect associated with the induced 
moment. We find that the induced magnetic moment 
strongly enhances the local quasiparticle density of states 
near the impurity. In addition, the blocking effect of the 
BiO surface layer between the STM tip and the probed 
CuO plane is essential to explain the observed resonance 
peak on the Zn site and the spatial pattern of the dif- 
ferential tunneling conductance around a Zn impurity in 
BSCCO. 
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FIG. 1. The local density of states on the near- 
est-neighboring sites of the unitary impurity A = (1,0) (a), 
B = (0,1) (b), and C = (-1,0) (c). Also shown with the 
dashed line the LDOS on the corner site of the supercell, 
which resembles the bulk DOS for a clean d-wave supercon- 
ductor. Here the A site is the position of the induced moment 
with the nonmagnetic unitary impurity located at (0,0)-site. 
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FIG. 2. The local density of states on the near- 
est-neighboring site of the nonmagnetic unitary impurity after 
the average over the four moment-position configurations. As 
a comparison, also shown with the dashed line the LDOS on 
the corresponding site for a single nonmagnetic unitary im- 
purity (i.e., without induced magnetic moment). 



FIG. 3. A two-dimensional bubble view of the spatial dis- 
tribution of the local density of states at the resonance energy 
E — 0.02 around the nonmagnetic unitary impurity with a 
magnetic moment. The bubble size on each Cu site scales the 
LDOS intensity. 
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